Decaying Vacuum Cosmology and its Scalar Field Description 
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We discuss the cosmological consequences of an interacting model in the dark sector in which the 
A component evolves as a truncated power series of the Hubble parameter. In order to constrain 
the free parameters of the model we carry out a joint statistical analysis involving observational 
data from current type la supernovae, recent estimates of the cosmic microwave background shift 
parameter and baryon acoustic oscillations measurements. Finally, we adopt a theoretical method 
to derive the coupled scalar field version for this non-equilibrium decaying vacuum accelerating 
cosmology. 

PACS numbers; 



I. INTRODUCTION 

In the last decade, the analysis and interpretation 
of observational data are strongly indicating that the 
present Universe is accelerating and that its geometry 
is spatially flat [ll-U- These observations gave rise to a 
resurgence of interest in a non-zero cosmological constant 
A that dominates the current composition of the Universe 
0. The main motivation behind this idea is that a very 
small value for A (~ 70% of the cosmic composition), not 
only would explain the observed late-time acceleration 
of the Universe (as indicated by type la supernovae ob- 
servations but would also reconcile the inflationary flat- 
ness prediction, corroborated by cosmic microwave back- 
ground data, (r^Totai — 1) with the current clustering 
estimates that point systematically to fi™ ~ 0.3 [l, @. 

Although flat models with a very small cosmological 
term are in good agreement with almost all sets of obser- 
vational data, the value of A inferred from observations 
(pa = A/SttG ~ lO-^^GeV) is between 50-120 orders of 
magnitude below estimates given by quantum field theory 
(pa ~ 10''"'^ GeV*). This extreme discrepancy originates 
an extreme fine-tuning problem, the so-called cosmolog- 
ical constant problem, and requires a complete cancella- 
tion from an unknown physical mechanism (Tj-fllj . 

Another problem with the cosmological term is that, 
although a very small (but non-zero) value for A could 
conceivably be explained by some unknown physical sym- 
metry being broken by a small amount, one should be 
able to explain not only why it is so small but also why 
it is exactly the right value that is just beginning to 
dominate the energy density of the Universe now. Since 
both components (dark matter and cosmological term) 
are usually assumed to be independent and, therefore, 
scale in different ways, this would require an unbeliev- 
able co-occurrence, usually referred to as the coincidence 



problem (CP). 

Many phenomenological models with variable cosmo- 
logical term (decaying vacuum) have been proposed in 
literature as an attempt to alleviate the cosmological con- 
stant problem |K)|-[13| and more recently the coincidence 
problem [13l - ll6| . In the context of the general relativity 
theory a cosmological term that varies in space or time 
requires a coupling with some other cosmic component, 
so that the total energy- momentum tensor is conserved. 
Thus, A{t) models provide either a process of a particle 
production or a time-varying mass of the dark matter 
particles increases [l7| . 

In this paper, we focus our attention on the theoretical 
and observational aspects of a A(i)CDM model in which 
the cosmological term evolves as a truncated power series 
of the Hubble parameter. The decaying vacuum compo- 
nent is assumed to be coupled only with the dark matter 
(interacting dark sector) thereby producing dark matter 
particles. The work is structured as follows. In Section II, 
we discuss the general aspects and solve the basic equa- 
tions of the decaying vacuum model. The observational 
constraints on parameters of the model are discussed in 
Section III. In section IV, by following an approach pro- 
posed by Maia and Lima |18|, we also derive a coupled 
scalar field version for this vacuum decay scenario. Fi- 
nally, the basic results are summarized in the conclusion 
Section V. 



II. THE A{t) MODEL 

Let us first consider a homogeneous, isotropic and 
spatially flat universe described by the Friedmann- 
Robertson- Walker line element. In such a background, 
the Einstein field equations are given by 



and 
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BttGp/ + A 3i?^ 



SnGpf - A ^ -2H - 3H^ 



(1) 



(2) 



where pf and pf are, respectively, the energy density and 
pressure of the cosmic fluid (radiation, baryons and dark 
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matter) and the dot denotes the derivative with respect 
to time. 

In this paper, we wih work with the A(t) function writ- 
ten as a power series of the Hubble parameter up to the 
second order, i.e.. 



A{H) ^ X + aH + 'SI3H^ 



(3) 



where A and a are constants with dimensions of and 
H, respectively, while /3 is a dimensionless constant. The 
term proportional to was proposed long ago by Car- 
valho, Lima and Waga [10^ based on dimensional argu- 
ments. Later on, it was also justified by using techniques 
from the renormalization group approach [l9| . The linear 
term was first discussed by Carneiro and collaborators 
[20j . The interest in the extended model represented by 
the above truncated power series was also suggested in 
the appendix C of a paper by Basilakos, Plionis and Sola 
2l|. As remarked before, here we discuss the observa- 
tional potentialities of this model, and, for completeness, 
we also derive the associated scalar field description. 

Now by combining the above equations and considering 
that pf = (7 — l)pf where 7 is the adiabatic index, we 
find 



2H + 37(1 - I3)H^ - ajH - 7A = 



(4) 



The integration of above equation leads the following H 
parameter 



H = - = , 
a 3(1 



1 



1 



7 e" 



1 



(5) 



where a = {'^/2)y^a^ — 12A(1 — (3) and a is the cosmic 
scale factor. In what follows, we will restrict our analysis 
to case in which A = 0. The practical reason for this 
choice is that the linear term in Eq. Q causes a tran- 
sition decelerating/accelerating phase. So that, we can 
obtain the following solution for the scale factor 



a{t) = C7(e^^*/2 - 1)2/37(1-/3) 



(6) 



where C is a integration constant. Combining Eqs. ([5]) 
and ([6]) one finds an expression for H{a), i.e.. 



H = 



1+ - 



3(l-/3)/2- 



(7) 



and we have used 7 = 1 in the last equation. By tak- 
ing Eqs. (H)), ([3]) and ([7]) at current time it is possi- 
ble to show that C = n^/{l - 13 - 0^)2/3(1-'^) and 
a = 3Ho{l — 13 — fim), where f2„j — SirGpmfi/iHQ. Now, 
by considering that a = (1 + z)^^, the Hubble parameter 
can be rewritten as 



H = H. 



1 - 



+ ^(l + ^)3(l-/3)/2 



1-/3 1-/3 



(8) 



It should be noticed that for /3 = 0, the above expression 
reduces to the dynamical A solution derived in Ref . [l^l . 



It is also straightforward to show from (|S]) that the 
deceleration parameter, defined as q{a) = —ad/a^, now 
takes the following form: 



1 



(9) 



Note that for z = 0, one finds qo = 1.5ilm — 1, and, 
therefore, the present value of the deceleration parame- 
ter is negative (see Table 1) and does not depend on the 
values of /3. It is easy to check that the transition red- 
shift (z(), defined to be the zero point of the deceleration 
parameter, is given by 



zt 



2(1 - /3 - ftra) 



(1 - 3/3)1], 



2/3(1-/3) 



(10) 



In this connection, we recall that some authors have 
recently suggested a picture based on theoretical and ob- 
servational evidences where the cosmic acceleration may 
have just peaked and started to slow down From ex- 
pression ([9|) for the deceleration parameter, we see that 
in the future (negative redshifts) the value of q(z) ap- 
proaches -1. 

For the sake of completeness, we derive the age-redshift 
relation for this class of A(t) models. We find 



t{z) = - In 

(7 



1 



1 ~ /3 - 



■(1 + ^) 



(11) 



Note also that the value of the present age of the Universe 
in this decaying vacuum model is given by 



2Hn 



3(l-/3-n„ 



■In 



1-/3 



(12) 



so that HqIq ^ 1 in agreement with the current observa- 
tions (see Table 1). 



III. OBSERVATIONAL ANALYSIS 

In this section we will discuss bounds on the free pa- 
rameters and /3. To this end we will use different 
observational sets of data, as described below. 

A. Constraints from Type la supernovae 

In this test we use one of the most recent SNe la 
compilation, the so-called Union 2.1 sample compiled in 
Ref. [i^l which includes 580 data points after selection 
cuts. 

The best fit to the set of parameters s = (ilm, /3) is 
found by using a statistics, i.e.. 



^ K(flf)(f«) 



(13) 
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where ^p{z\s) — Slogd^ + 25 is the predicted distance 
modulus for a supernova at redshift (z), cIl is the lumi- 
nosity distance, iJ,l,{z\s) is the extinction corrected dis- 
tance modulus for a given SNe la at Zi and cTi is the 
uncertainty in the individual distance moduli. 

B. Constraints from CMB/BAO ratio 

Additionally, we also use measurements derived 
from the product of the CMB acoustic scale £a = 
7rdyi(2:*)/rs(z*) and from the ratio of the sound hori- 
zon scale at the drag epoch to the BAO dilation scale, 
rs{zd) / Dv{zBAo), where dA{z*) is the comoving angular- 
diameter distance to recombination (z* = 1089) and 
rs(z*) is the comoving sound horizon at photon decou- 
pling. In the above expressions, Zd — 1020 is the red- 
shift of the drag epoch (at which the acoustic oscilla- 
tions are frozen in) and the dilation scale, Dy, is given 
by Dv{z) — [zr'^{z)/H{z)]^^^. By combining the ratio 
r,{zd = 1020)/r/(z* = 1090) = 1.044 ± 0.019 0, 
with the measurements of Vg (zd) / Dy (zbao ) at zbao = 
0.20, 0.35 and 0.6., one finds [Ull^ 

/o.2o = dA{z,)lDv{Q.2) = 18.32 ± 0.59 , 
/o.35 = rfA(^*)/^y(0.35) = 10.55 ± 0.35 , 
/o.6 = dA{z.,)/Dv{Qm) = 6.65 ± 0.32 . 




FIG. 1: The results of our statistical analysis. The constraints 
from SNe la and CMB/BAO ratio data are shown by shaded 
contours for Scr confidence level, while the contours from com- 
bining observational data are represented with black lines for 
(68.3% CL, 95.4% and 99.73%) confidence intervals. 



C. Combining observational sets data 



Let us now discuss the constraints given by the above 
different sets of data. In our statistical analysis, we mini- 
mize the function x\ = Xsn + Xcmb/bao' where xln ^nd 
^CMB/BAO corresponding to the SNe la and CMB/BAO 
functions, respectively. 

In Figure 1, we display the space parameter (f2m,/3). 
The supernova data set fitted with Union 2.1 are shown 
by the bluish contours while the orange contours are 
derived from CMB/BAO ratio data alone. The re- 
sults from our joint analysis are shown by dashed lines. 
By marginalizing on the nuisance parameter h (Hq — 
imhkms-Hdpc-^) we find /3 = -0.185t°;°^^l° °^™f5 
and - 0.325lg;g^^ig;g^^lg;i;^;5 at 68.3%, 95.4% and 
99.7% of confidence level, respectively, with x^j„ = 564 
and V = 581 degrees of freedom. The reduced Xr — ^■^'^ 
where (x^ = Xmin/'^)^ thereby showing that the model 
provides a very good fit to these data. We also note 
that, due to the process of matter production resulting 
from the vacuum decay, the current observational bounds 
on il„i provide a value slightly higher than the ACDM 
model. 

By using the best-fit values of ilm and /3 it follows 
that qo = -0.513;!;°^i^, zt = 0.991°;^^! and toHo = 



TABLE I: The results from joint statistical analysis 



Parameter best fit error bars tj^" 





0.325 


+0.035 
-0.035 




-0.185 


+0.045 
-0.055 


a/Ho 


2.58 


+0.24 
-0.27 




-0.513 


+0.053 
-0.053 


Zt 


0.99 


+0.19 
-0.20 


to Ho 


1.003 


+0.049 
-0.051 



+0.049 
-0.051- 



1.003 
shown in Table 1. 



The main results of our joint analysis are 



IV. SCALAR FIELD DESCRIPTION 

A(t) cosmologies constitute a possible way to ad- 
dress the cosmological constant problem. However, in 
the absence of a natural guidance from fundamental 
physics, one needs to specify a phenomenological time- 
dependence for A in order to establish a definite model 
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and study their observational and theoretical implica- 
tions. A possible way to seek for physically motivated 
models is to represent them through a field theoretical 
language, the easiest way being through scalar fields. In 
this section, we discuss a possible derivation of a time- 
varying A models from fundamental physics based on the 
approach suggested in Rcf. [18)] (see also (2^] for a general 
description with arbitrary curvature). 

Following standard lines [l^ , let us define the function 



7* 



2H 
MP 



= 1 - 



A 
MP 



(14) 



Each K{t) model is specified through the 7^ time- 
dependent parameter. Thus, for the vacuum decay sce- 
nario presented here, the function 7^ takes the following 
form: 



7* = 



3(l-/3)i?2 -cfH 
3H^ 



(15) 



As an intermediate step, it is convenient to replace the 
vacuum energy density and pressure as given by Eqs. ^ 
and ^ by the corresponding scalar field expressions, i.e.. 



A/SvrG ^ P0 - + V{cf>) , 



and 



- A/SttG 



^2/2 - V{cj,) 



(16) 



(17) 



where and are, respectively, the energy density and 
pressure associated to the coupled scalar field </> whose 
potential is F(0). 

Now, by defining a parameter x = 0^/ (0^ + pdm) with 
< X < 1, we can manipulate Eqs. ([1]), ([2]), (fTHl) and 
([TT]) to separate the scalar field contributions (see for 
more details), i.e.. 



and 



y(</>) 



3IP 
8^ 



3IP 
8^ 



(18) 



(19) 



Note that the above equations link directly the field 
and its potential with the related quantities of the dy- 
namical A(t) case. From Eq. (fT8)) . one can show that the 
field (j) in terms of the Hubble parameter reads 



± 



SttG 



Ho 



H 



Now, inserting Eq. ([T5|) into Eq. (|20p and integrating 
it, one obtains 



= -A-Mn 



(21) 




FIG. 2: Behavior of the scalar field potential in units of the 
present critical energy density U{<j)) = V{(j))/pc,o- The dotted, 
solid and dashed curves correspond to the lower, best fit and 
upper limits of the Qm and /? parameters at {la C.L.). See 
discussion in the main text and Table 1. 



where A = y/3TTG{l ~ l3)/2x and A = cr/3(l - /3). By 
using Eq. (|19p one finds 



647rG 



[(2-x)a + i?/(<^)]/(^) , (22) 



where 
/(0) 



A y 



g -2A(0-0o) ^ [ _ ) g2A(</,-</,o) _^ 2 



D2 / 



A 



, (23) 



with B = 3D^[2/3+{l-P)x]/4: and D = y/Ho - A+y/Ho. 
Note that by assuming cr = A — 0, and the po- 
tential takes the form of a simple exponential function 
thereby reducing to the potential initially found in Ref. 
fisj. It is worth mentioning that double exponential po- 
tentials of the type (|22p have been considered in the lit- 
erature as viable examples of quin tessence scenarios (see, 
e.g, [28|]). As discussed in Ref. [23], a scalar field potential 
given by the sum of two exponential terms is also moti- 
vated by dimensional reduction in M-theory with inter- 
esting implications for the late-time accelerating behavior 
of the cosmic expansion. 

In Figure 2, we display the evolution of the scalar 
field potential (in units of the present critical energy 
density [U{(j)) = V {(j)) / pco]) for some selected values of 
the free parameters, namely: ftm = 0.29, 0.325, 0.36 and 
13 = -0.24, -0.185, -0.14. This choice correspond to the 
lower, best fit and upper limits on rim and /3 parame- 
ters at (1(7 C.L.). Note that all curves exhibit the same 
behavior differing only by a small shift. 
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V. CONCLUSIONS 

In this paper we have investigated theoretical and ob- 
servational features of a vacuum decay model in which 
the cosmological term is written as a power series of the 
Hubble parameter up to the second order. We have per- 
formed a statistical analysis involving the lates obser- 
vational measurements of SNe la, BAO peak and CMB 
shift parameter and found stronger constraints on the 
parametric space rim — P- By using the best-fit values 
for = 0.325lJ];!]|^ and (3 = -0.185t|]:[!^^ we have also 
derived the current values of go j zt and to , that are go = 
-0.513lg °i, zt = 0.99tS]ii^ and hH^ = 1.003t°:S]^?, re- 
spectively. These values are in excellent agreement with 
the predictions of the standard model (see Table 1). 

By following a theoretical method developed in Ref. 
[isj we have also derived the scalar field description for 
this cosmology and shown that the present A(t)CDM 



model is identified with a coupled quintessence field with 
the potential of Eq. (P^ . Since this A(t)CDM scenario is 
quite general (having many of the previous proposals as 
a particular case), we argue that a coupled quintessence 
field model whose potential is given by Eq. ([^ is dy- 
namically equivalent to a large number of decaying vac- 
uum scenarios previously discussed in the literature. 
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